Comparative Analysis of Pion, Kaon and Proton Spectra Produced at PHENIX 
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The shapes of invariant differential cross section for identified w ,K , p and p production as 
function of transverse momentum measured in pp collisions by the PHENIX detector are analyzed. 
Simultaneous fit of these data to the sum of exponential and power-law terms show significant 
difference in the exponential term contributions. This effect qualitatively explains the observed 
shape of the experimental K/tt and p/ir yield ratios measured as function of transverse momentum 
of produced hadrons. A picture with two types of mechanisms for hadron production is given. 
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There is a large volume of experimental data available 
on charged particle and identified hadrons production in 
high energy particle collisions collected starting from the 
very first experiments performed on the Intersecting Stor- 
age Rings (ISR) at CERN and till contemporary high 
statistics measurements carried out at RHIC and LHC. 
In most available publications the measured experimen- 
tal spectra of produced hadrons become a subject for a 
phenomenological description or for a comparison with 
selected models in each experiment separately. However, 
a comparative simultaneous analysis of the whole avail- 
able data volume could provide a new powerful lever arm 
to disclose a common underlying dynamics in work in 
the hadron production in high energy particle collisions. 
Such systematic analysis have been described in two our 
recent papers [TJ [2] . 

In pQ it was demonstrated that spectra of charged 
hadrons produced in collisions of baryons require an exis- 
tence of a sizable fraction of an exponential (Boltzmann- 
like) statistical ensemble of charged hadrons on the top of 
a power-law (pQCD inspired spectrum shape) functional 
term. According to [T] the overall generic charged hadron 
spectrum as function of the produced hadron transverse 
momentum (Pt) is given by a sum of the exponential and 
power-law terms 
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where Exkin = V Pj- + M 2 — M with M equal to the 
produced hadron mass. A e ,A, T e ,T, n are the free pa- 
rameters to be determined by fit to the data. 

A typical charged particle spectrum as function of 
transverse energy, fitted to the function (nj) is shown in 
Fig [I] In addition the contributions off the exponen- 
tial and power law terms to the spectrum are shown in 
Fig [T] separately. It is observed, that the exponential 
Boltzman-like term dominates the charge particle spec- 
trum at low Pt values. Since the absolute majority of 
produced charged hadrons are charged tt it was concluded 
that the observed charged hadron spectrum shape is also 
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FIG. 1. Charged particle spectrum [3] fitted to the func- 
tion ([l]): the red (dashed) curve shows the exponential term 
and the green (solid) one stands for the power law term. 



characteristic for charged tt. This conclusion was sup- 
ported by the studies of the identified charged tt produc- 
tion spectra pQ. 

Though the data on the identified hadron production 
are sparse and available in the very limited range of the 
hadron transverse momentum only it is still interesting 
to analyse the spectra of K and p in a way proposed 
in pp. This analysis was performed in 2\. It turned 
out that, contrary to the pions, the spectra of K and 
p leave a relatively little room for the exponential term 
contribution and within the experimental errors could be 
described by the power law term only. 

In the present paper an attempt was made to make a 
simultaneous comparative analysis of the it,K and p spec- 
tra produced at the same collision energy and under the 
same experimental conditions. For this comparison the 
data from the detector PHENIX collected in pp collision 
run at y^s) = 200GeV were used [J. 

As the first step one could note a peculiar behavior of 
the ratios of differential cross sections K/tt and p/ir mea- 
sured in the experiment. These ratios plotted as function 
of Pt of produced hadrons are shown in Fig [2^,, b. In 
both cases the ratio reaches a plateau above Pt ~ 2GeV 
and drops down for low Pt values. Within the frame- 
work of the proposed approach based on the formula (fTl) 
the observation of this plateau suggests that parameter 
n of the power law term in M is likely to be similar both 
for K, p and it. In the QCD model this is correct at 



Type of 
produced hadron 


Exponential 
contribution [%] 


Power law 
Terms 


Exponential 
Terms 


TV 


75 ±2 


- 


- 


K 


28 ±7 


0.29 ±0.03 


0.04 ±0.003 


V 


18 ±8 


0.25 ±0.02 


0.02 ±0.001 



TABLE I. Contribution of the exponential (Boltzmann-like) 
term to the charge particle spectra and K/ir and p/ir val- 
ues calculated separately for power-law and exponential terms 
of |l| obtained from the simultaneous fit of ir,K and p [J. 



high Pt values since the produced hadron distributions 
are largely driven by the gluon momentum distribution 
in the colliding particles. If so, the identified ir,K and 
p spectra are fitted simultaneously to the function 
with a constraint to have the same value for n parame- 
ter in the power law term for all three spectra together. 
As the measured charged kaon spectrum is restricted to 
the low- Pt values it was extended using the available K® 
data. The result of this fit procedure is shown in Fig [3] 
and table [I] It is important to note that the values for 
the parameter T obtained from such fit procedure turn 
out to be practically the same for all types of produced 
hadrons. Such surprising result was not obvious a priori. 
As it was expected the Boltzmann-like term for 7r produc- 
tion dominates, while it gives much lower contribution in 
K spectrum and it is practically close to zero for the p 
spectrum. The relative amount of the exponential term 
contribution to the hadron spectra estimated from the si- 
multaneous fit of TTjK and p differential cross sections to 
the function |l]) are given in Table 1. It is important to 
note that the observed significant difference in the expo- 
nential contributions to the hadron spectra implies the 
difference in the hadron production mechanisms rather 
than an artifact of the fit procedure. 

Once arrived to this conclusion, one could explain the 
peculiar behavior of the K/n and p/n ratios (Fig [5]) as 
function of Pt of produced particles. Figure [4] shows that 
the ratios of differential cross sections K/ir and p/ir mea- 



sured in the experiment and those obtained from the fit 
are in a good agreement with each other. The existence of 
large fraction of low-P T tt represented by the Boltzmann 
exponential statistical distribution suppresses these ra- 
tios at low values of Pt- This situation is somewhat 
similar to that found in |T] for tt production in pp and 
heavy ion collisions taken at the same collision energy 
per nucleon. The ratio of the correspondent differential 
cross sections as function of Pt significantly differs from 
unit and is traditionally interpreted as a signature of the 
nuclear absorption effect. It has been shown in [T] that 
the shape of this ratio is formally defined by a difference 
of partial contribution of the exponential term to the 
sum (IT]) for minimum bias Au — Au and pp interactions. 
Within the framework of the proposed approach the 
hadron production mechanism could be represented in a 
way of a sum of two different contributions: 

1) release of quasi-thermalized hadrons (mostly n, and 
much less probable K and p), described by the exponen- 
tial Boltzmann spectrum and 

2) pQCD like production of hadrons described by the 
power law statistical distribution. The parameters of this 
distribution do not depend on the type of hadron pro- 
duced in pp collisions. 

It is worthwhile to note that as was found in [1] a siz- 
able contribution of the thermalized hadrons shows up in 
the baryonic collisions only. Whilst the interactions with 
the high energy photons involved as colliding particles do 
not require an extra exponential statistical distribution 
to describe the produced hadronic spectra. 

In conclusion, it is found that the 7r production in pp 
collisions is dominated by a release of quasi-thermalized 
particles, while the spectra of heavier K, p and p are 
dominated by pQCD like production mechanism, leaving 
a relatively small room for thermalized particle produc- 
tion. For it, K and p production the parameters of the 
power law contribution take practically the same (within 
the errors) values which depend on the global conditions 
like a collision energy, or a type of colliding particles 
rather than a type of produced hadron. 
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FIG. 2. Ratios of the measured differential cross sections K/tt (a) and p/iv (b) measured at PHENIX [4]. For (a): K - blue 
points K® - green points. For (b): p/iv~ - blue points, p/-7r - green points. 
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FIG. 3. n, K and p spectra U fitted simultaneously: the red line shows the exponential term and the green one stands for 
the power law. The exponential term dominates for 7r spectra only. 
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FIG. 4. K/n and p/n ratios for experimental data [4] (points) shown over the corresponding ratios obtained from the results 
of the simultaneous fit (line). The error band is calculated using the parameter uncertainties obtained from the fit. 



